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ABSTRACT 


The application of numerical optimization techniques to 
the design of drive shafts is demonstrated. The analysis 
investigated the effects of a small mass imbalance in con- 
junction with the rotation of a shaft with synchronous 


whirl. 
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T OPP PPODUCTION 


Optimization and composite materials are two fields 
that have seen considerable developments in recent vears. 
Advancements in these technologies have provided the design 
engineer with new dimensions, the depths of which are yet 
to be explored. 

The major material characteristics that are advanta- 
geously cultivated in composites are the high strength and 
stiffness to weight ratios. “he weight reduction is so 
significant that the aircraft and aerospace industries have 
concentrated a major portion of their reserach efforts in 
this field. Investigations in improving the resistance of 
composites against environmental factors such as temperatures, 
corrosion, and wear are continuing and findings affirm the 
contention that the future of this field is bright. However, 
other industries are slow in utilizing these materials in 
basic machinery elements, mainly due to a major drawback-- 
the high cost of manufacturing and processing. It is in 
this vein that the importance of an optimized design becomes 
apparent.  Advancements in fabrication and processing of 
composite materials coupled with the capability to optimize 
designs may well lead composites into a more competitive 


market. 





! 
! 
! 


— — 


With the ultimate goal of optimizing the design of a 
drive shaft made of composite materials, a research program 
was initiated. This report is the first in this investiga- 
tive effort. 

Before the intricacies introduced Бу the directional 
material properties of composites are included in the 
analysis, it was deemed necessary to probe the characteris- 
tics of the design of an isotropic material. This will 
enable one to foresee the trends and the adaptive measures 
necessary when the composite material properties are con- 
Sidered. The application of numerical optimization tech- 
nigues to the design of an isotropic drive shaft is demon- 
strated here. The effects of a small mass imbalance on the 
bending moment and deflection of a shaft with synchronous 
whirl are investigated. Mainly, pose. the objective is 
to develop an analysis procedure that may be adapted to the 
peculiarities of composites and ultimately be used in op- 
timizing a composite drive shaft design. 

Section 2 defines the scope and limitations of the 
analysis and optimization. 

Section 3 discusses the analysis itself, the equations 
and the modes of failure that were considered. Examples of 
the analysis are shown. 

Section 4 briefly describes the optimization program 


[1,2] and presents the results of optimizing a number of 





shaft designs. “wo variable function space diagrams are 
provided to illustrate the design space for different load 
Conditions. 

Section 5 contains recommendations for future investi- 
gations. 

Appendix A contains a description of the FORTRAN program 
used in the analysis and the program text. Appendix B shows 
the derivation of equations for deflection, moment and shear 


force. 
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πι SCORE SANDS LIMITATIONS 


The analysis applies to shafts made of isotropic ma- 
terials with simply supported or fixed-fixed end conditions. 
Basic shaft design formulas are used except that the de- 
flection and moment equations are derived to include the 
effects of a small mass imbalance. The shaft is assumed 
to rotate with synchronous whirl. 

To insure a continuous design field, the optimization 
1s constrained to subcritical speeds only. Supercritical 
speeds introduce feasible design regions that are disjoint 
from the primary design space. In this investigation, these 
secondary areas are considered infeasible. 

Since the optimization tends to a large radius-small 
thickness design, radial stresses are neglected. However, 
buckling of thin cylinders due to torsion and compression 
is considered. The equations used for these failure criteria 
do not include dynamic effects. A constant torque and a 
Steady, uniform axial force are the loads considered in 
the formulas [3]. 

The analysis incorporates the capability to design a 
shaft that may be used in two or more loading conditions. 
For example, a shaft may be designed to transmit 50 HP at 
300 RPM with an axial load of 2000 lbs. as well as to transmit 
150 HP at 3000 RPM with a 1000 1b. axial load. 


11. 





III. ANALYSIS 


The problem is a shaft required to transmit a specified 
horsepower at a given speed. There may be an axial load, F 
or an internal presence, P. (See Figure 1.) Basic equations 
are used throughout the analysis and fundamentals such as 
formulas for area, volume, moment of inertia, polar moment 
of inertia and radius of gyration are not repeated here. 
However, certain stress calculations will be emphasized to 
illustrate the flow of analysis. 

Appendix B shows the derivations of the equations for 
bending moment, deflection and shear for pinned-pinned and 
clamped-clamped end conditions. These derivations have 
taken into consideration the effects of a small mass imbalance 
in conjunction with the rotation of a shaft with synchronous 
whirl. 

The inch-pound-second system of units (IPS) is used in 
the computer program analysis. However, bv changing the 
equation for torque (line 50 in the FORTRAN text, Appendix A) 
to reflect the International System of Units, the program 


may be used effectively with consistent SI units. 


η 








Figure 1. Applied Loads: Axial load (F), Internal pressure (P), 
Torque (T), Mass imbalance (e). 


a.  Pinned-pinned end condition 


ም 
> A m 
SA ү 


b. Clamped-clamped end condition 


= l- cos —— 


Šo [ ли) 
2 ዴ 


Figure 2. Assumed deflected shape of the shaft. 
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A. SIMPLY SUPPORTED BEAM 
The deflected shape of the beam is assumed as: 
6 = on sin > (see Fig. 2). The solution of the govern- 


ing equations equations is given in Appendix B. The deflection 


is maximum at x = 2/2 
E |. a Cae 
max 384 (2 i 


4 
EI-K, ( = Е 


where: K. = DAW? 


o = mass per unit volume 

A = cross sectional area 

w = shaft speed is radians per second 
K = pAg = weight per unit length 


e = eccentricity of mass with respect to the axis 
of rotation 


2 = shaft length 
E = Young's Modulus 
I = Moment of Inertia 
F = Axial Load 
Certain characteristics of this equation are of particular 
interest. If there is no rotation nor axial load, the formula 
reduces to the classical equation for maximum beam deflection. 


ሙ፳ሓዴ 
max = (3.1b) 
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If there is no axial load, the equation becomes unstable 
(deflection goes to infinity) when EI = (5) Or Esel zung 


for the critical speed, 


[ο τσ 


which defines the first fundamental frequency of a simply 


supported shatt. 


If there is no rotation, instability arises when EI = Β΄ 
er 
TEI ። 
F. E (3:1a) 


recognizable as Euler's column buckling criteria. 
With both rotation and axial load, the critical speed 


becomes 


(3.1e) 





To avoid the instability regions, the denominator of 
equation (3.1a) must be given a lower bound greater than zero 
(point A in Fig. 3). Additionally, a maximum deflection has 
to be imposed to prevent computer overflow (point B in Fig. 3). 
These limits have no relevance with the deflection constraint 
used in the design optimization process (point C in Fig. 3), 
but are used only to prevent numerical ill-conditioning on 


the computer. 
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Beflection 
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Figure 3. Limits on the deflection. The solid curve repre- 
| sents the deflection values used in the analysis. 
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The bending moment corresponding to the deflection calcu- 


lated by equation (3.1a) at x - 2/. 15 
K д2 
е 


Е е R 2 2 
DL CU Ша а (3, 1፻) 





When the weight tends to cancel the deflection, the corre- 


sponding moment at the same cross section becomes 














K, 2° 
min š М дах NC . ᾱ- 
The shear force is maxlmum at x = 0. 
Io К.е? КЕ 
anm 1. 2 
ከ =- 5 (So 
B.  CLAMPED-CLAMPED BEAM 
The deflected shape of the beam is assumed as: 
= 21% i | | 
Ó = > [ 1 - cos — Е (See Fig. 2b.) The deflection is 
maximum at x = çs 
2224 Bed > (3.2a) 
384 EI - 0.746 K,2 - 96 (=) Е 


As in the simply supported beam, similar observations 
regarding special loading conditions are apparent. If there 
is neither rotation nor axial load, the formula for maximum 


beam deflection results 


4 
ጅጋ“ 


max 384 EI 


О» 


(9. 29) 


E7 





With no axial load, instability occurs when 


ИСЕ - 0 746 ρλω΄ or 


= 2275609 — (5. ο) 
ዕልዷ 


which is the first fundamental frequency for a fixed-fixed 
beam baseä on the assumeä deflection used here. 


When there is no rotation, instability arises when 


2 
384 EI = 96 E ОГ 
Т 
2 
F = ευ (3.2 
С ¿2 


Buler's column buckling criteria. 
With rotation and axial load, the instability occurs at 


a smaller speed than that given by equation (3.2c), 






EI - 0.02533 r^ 


PAL” 






(3728) 


The bending moment at x - P corresponding to the 


deflection given by equation (3.2a) is 


2 


К. 7, К. 
El 52 36 2 
DEL | = + == + ELMO PS + = (3,2፻) 





A lesser bending moment occurs at the same cross section 


when the weight tends to cancel the deflection. 
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K, 
ti) 2:0 12. eden 
The maximum shear occurs at x = 0. 


κε ον К.е? K. 
EI j 2 
s. = 1 T 5 + = (3.2h) 








The preceding calculations are performed in Subroutine 


BEND, given in Appendix A. 


С. STRESS ANALYSIS 

With the deflection, shear and bending moments known at 
the critical cross sections, the analysis proceeds with the 
calculation of the stresses at the critical stress elements. 
At midspan, the critical stress elements are on the outer 
surface of the shaft where the axial stress (tensile or 
compressive) is maximum (Points A and B, Fig. 4). Since the 
shaft is in synchronous whirl, the same stress element, A, 
will generally tend to be in tension throughout the rotation, 
with the magnitude of the axial stress varying as the weight 
effect becomes alternatingly additive and subtractive. 
Element A is in the direction of the mass eccentricity from 
the center of the shaft while element B is on the outer 
surface in the opposite direction. 

Most of the following equations are found in basic design 
and applied mechanics books. References will be noted when 


necessary. 
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Axlal stress due to bending moment (tension/compression) 


መ J ο 

б. (тах) E (3.38) 
Puna To 

o (min) = = 


Axial stress due to axial load (compression) 


= 
2 (3. 3b) 


The axial stress due to internal pressure for thin-wall 


cylinders is: 
G = On (σοι ος] 


Superposition with careful regard to sign conventions 
provides the maximum and minimum tensile and compressive 
axial stresses. 

Hoop stress is primarily due to internal pressure and 


rotational inertia [4] 


μα. 
EDS = E (3730) 
2.2 r. 
PSS Z о asi] (3.3e) 
a ee ТЕЛЕ 
Н g 4 О 


where Y = weight per unit volume 
w = rotation in radians per second 
v = Poisson's ratio 
r. = inside radius 
r = outside radius 
g = acceleration due to gravity 


20 





The shear stress at x - AP is caused only by torsion 


(shear force is zero). 


о О 
TT = 77 USE) 
where Тт = shear stress due to torsion 
T = torque 
J = polar moment of inertia 


At x = 0, the critical element is where the shear stress 
due to the shear force along the cross section is additive 
with respect to the shear stress due to torsion (element C, 


Fig. 3). The shear stress due to um for thin wall cylinders 


ax 
Bl is 


D. max (3.39) 


Superposition of equations (3.3f) and (3.39) gives the total 
shear stress at the critical stress element. 

The axial stress at element C due to bending moment 38 
Zero. The axial stress at this point is, therefore, caused 
only by the axial load and the internal pressure, equations 
1) апа (3.3c). 

The hoop stress is calculated by equations (3.3d) and 


(3.3e). 


21 





stress element А 


stress element D 


a) At midspan 


Stress element C 





Оң (hoop stress) 






T (shear stress) 
D 


ር (axial stress) 


c) Typical stress element. Radial 
stresses are neglected. 


Figure 4. Critical Stress Elements 
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D. MODES OF FAILURE 
Eight modes of failure are considered in this investiga- 


tion. The first two evaluate the static strengths at the 


critical stress elements using Von Mises criteria [5]. At 
midspan, 
2 2 2 
σι 7 [Og ση or + o. + 3 Тт (3.4a) 
At x = 0 
Din 2 2 2 
m = Sy Sy С. + en + 3 Т ұрыну (3.4b) 


Failure is assumed if the yield strength of the material is 
equal to or less than these Von Mises stresses multiplied by 
a safety factor. 

A limit on the deflection is imposed. A common practice 
is to specify a maximum deflection in inch per foot of shaft 


length. For example, as used in this report, 


SH = 0.005 inch/foot length 


Since the optimization may lead to a thin cylinder design, 
buckling due to torsion and compression are considered. The 
critical torque and critical compressive stress for buckling 


ofthin cylinders [3] are; 


Q7 li E. . E (3.4c) 
S) (Esch eth ` 
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с _ _ Et _ (3.44) 


> r 3 (1:2) 


Experimental data indicate that actual failure of the cylinders 
usually occurs below fifty percent of the values calculated 
from these equations. To accommodate this discrepancy and 

the fact that the dynamic effects are not considered in the 
derivation of these equations, comparatively higher safety 
factors must be used for these failure modes. 

The shaft is also checked for column buckling using Johnson's 
or Euler's equations [5] as applicable. The actual slenderness 
ratio is compared with the slenderness ratio where transition 
occurs to determine which formulas apply. The calculated 
critical load, using the appropriate equations for the end 
condition, is compared to the actual load, again incorporating 
a safety factor. Equations (3.1d) or (3.2d) are used if the 
slenderness ratio falls in the Euler range. In the Johnson 


region, the following equations [5] apply. 


E 
e 


A 


፡ rn (ደነፆ (3.4e) 
= S b (к) 


where S = yield strength in compression 


(=) =s ]enderness ratio 
| 12 1. 
21 nE 
E - Young's Modulus 
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n = 1 for simple supports 
4 for fixed-fixed ends 
The effects of the fluctuating axial stresses are evaluated 
for fatigue failure using the Goodman diagram [5]. Von Mises 
stresses for the mean and alternating stresses are calculated 


as follows: 


Imax T min 
ο ሸሙ ጋ =s (541) 
max - “min 
С μη (3.4g) 
1 2 2 
On = fOr ση On + ση - 3 Тт (3.4Һ) 
om = m (3.41) 


The Goodman diagram is drawn (Fig. 5) and the safety factor 


is determined analytically. 





Figure 5. Goodman Diagram 
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SF = r = — (3.43) 
G t ! 
m 2 σ 
m a Se 
ο ο] 
σ t 
m 
where Se = fully corrected endurance limit 
Su = ultimate strength 


This safety factor is required to be greater than a specified 
value. 

Lastly, failure is assumed if a specified percentage of 
the critical speed as calculated by equations (3.le) or (3.2e) 


is less than the shaft speed. 


E. ANALYSIS EXAMPLES 

The examples are chosen to illustrate several aspects of 
the shaft design problem. The first is a design which vio- 
lates almost all of the constraints while the second example 
is one within the feasible region of the design field. 

The material used in all the examples in this report is 
UN G10150 HR steel, the properties of which are listed in 


Table I [5]. 
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TABLE I 
MATERIAL PROPERTIES 


UNST G10150 HR STEEL 


Young's Modulus 30.0E 06 psi 
Shear Modulus ΤΕ 06 psi 
Poisson's Ratio 005292 3 
Weight Per Unit Volume 0.282 lbs/in 
Yield Strength PSOE 03 psi 
Ultimate Strength 50.0E 03 psi 
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1. Analysis Example 1 
A simply supported shaft with 0.1 inch eccentricity 

is to transmit 150 HP at 3000 RPM (3150 in-lb torque), 
with an axial load of 2000 lbs. For a shaft length of 
120 inches, inside radius of 1.0 inch and thickness of 
1.0 inch, the design violates the static and dynamic 
strength, deflection, buckling due to compression, and 
the speed constraints. The analysis of this design is given 
below where an asterisk (*) denotes violated constraints on 


the design. 


Volume 0.11310Е 04 іп? 
Weight 0.31393E 03) lbs 
Moment of Inertia ` 0.11781E 02 in? 
Polar Moment of Inertia 0.23562E 02 in^ 
*Critical Sreed 0.14767E 04 RPM 
Critical Axial Load 0.18764E 06 155 
арса Bucklrng"wOrque 0.47469E 08 in lb 
*Critical Buckling Stress 0.905493 07 psi 
*Deflection 0.24000E 02 in 
Bending Moment l 0.36420E 08 in 1b 
Shear Force ШЕТПЕ Ос lb 
Axial Stress at À 0.95779E 07 psi 
Axial Stress at B -0.95783E 07 psi 
Hcop Stress (" 501057 03 psi 
Torsional Shear Stress 07267383E 03 psi 
Total Shear Stress 0.13337E 06 psi 
Mean Axial Stress 0957 UE 07 psi 
Alternating Axial Stress 0.81250E 03 psi 
*Von Mises at Midspan ο 75 5FE 07 psi 
*Von Mises at x = 0 0.23100E 06 psi 
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The calculated deflection is the maximum value 
imposed to avoid numerical ill-conditioning (point B, Fig. 3). 
Since the deflection affects most of the relevant parameters 
in the analysis, the calculated stresses are unreliable and 
not valid. These values, however, tend to direct the design 
towards the feasible region. 

2. Analysis Example 2 

This shaft is similar to the first example in design 
conditions and loading. The only difference is that the 
dimensions are such that no constraints are violated, but 
two constraints are active or critical. The shaft is 
120 inches long with an inside radius of 6.39 inches and 
thickness of 0.0369 inch. This, in fact, is the optimized 
design for the given loading as calculated by COPES-CONMIN 
[1,2]. The results of the analysis follows, where a double 


asterisk (**) denotes values at their design limit. 


Volume 0.17866E 03 in” 
Weight 0.50382E 02 lbs 
Moment of Inertia 0.30629E 02 in! 
Polar Moment of Inertia 0.61258E 02 in! 
**Critical Speed 0.60061E 04 RPM 
Critical Axial Load ΤΗΕ lbs 
**Critical Buckling Torque 0.31490E 05 in 1b 
Critical Buckling Stress 0.10400E 06 psi 
Deflection 0.43813E-01 in 
Bending Moment 0120268905 in lb 
Shear Force 0.84895E 03 lbs 
Axial Stress at A 0.22326E 04 psi 
Axial Stress at B -0.49193E 04 psi 
Hoop Stress 0.29778Е 04 psi 
Torsional Shear Stress О 50795 02 psi 
Total Shear Stress 0.14712E 04 psi 
Mean Axial Stress 0.20739E 04 psi 
Alternating Axial Stress 0.15872E 03 psi 
Von Mises at Midspan 0.69313E 04 psi 
Von Mises at x = 0 0.39193E 04 psi 
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In each of the examples, the analysis proceeds as 
follows. All shaft characteristics like cross sectional 
area, volume, weight, moment of inertia and polar moment 
of inertia are first calculated. Parameters that depend 
only on geometry such as critical torque and compressive 
stress for buckling of thin cylinders, and critical speed 
are then determined before Subroutine BEND, where the de- 
flection, bending moments and shear force at the critical 
Cross sections are computed, is called. 

At midspan, the absolute values of the axial 
stresses at elements A and B (see Fig. 4a) are compared 
to determine which element is critical. In both examples, 
element B is considered critical for strength evaluation. 
At x 2 0, the critical element, C, is where the torsional 
shear stress and the shear stress due to the shear force 
along the cross section are additive. The Von Mises 
stresses for each of the critical elements are then calcu- 
lated using equations (3.4a) and (3.4b). The rest of the 
analysis is a straight forward application of the equations 


for the critical parameters and the constraints. 
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ШУО КОРДИМ НУ ШОН, RESULTS AND CONCLUSIONS 


A. OPTIMIZATION AND RESULTS 

There are a number of optimization programs available, 
each one using different techniques in locating the desired 
optimum design.  COPES-CONMIN [1,2] is a versatile program 
that may be used for sensitivity analysis and two variable 
function space study as well as an optimization tool. 

Provided with a user supplied analysis program (Subroutine 
ANALIZ) where the objective function, constraints and other 
relevant parameters are calculated, it determines a usable 

and feasible sector from which a search direction is chosen. 
This choice of search direction is made by using such informa- 
tion as the gradients of the objective function, gradients of 
active and violated constraints and the degree of push-off 
these constraints violations generate. This iterative process 
of minimizing/maximizing the objective function by changing 
the design variables is terminated when no further improvement 
can be made. 

In this exercise, the objective function to be minimized 
is the volume, with the radius and thickness as design 
variables. Several design examples are shown to demonstrate 
the effects of the different parameters on the optimized 


design and on the design space itself. Seven examples are 
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presented corresponding to the loading conditions as listed 
in Table II. Two variable function space diagrams are 
generated to illustrate these effects. The numbered curves 
in the graphs represent the corresponding constraints shown 
in Table III. The shaded side of a curve represents the un- 
feasible area of the design space with respect to the cor- 
responding constraint. The dashed curves represent contours 
of the objective function, the numbers indicating the 


volume in cubic inches. 


TABLE II 


OPTIMIZATION EXAMPLES 


ES Horsepower RPM Pressure Axial Load 
(psi) (lbs) 
0 0 





*#Example No. 7 is a design for two load conditions. 
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Number 


ο ፌ] ርነ ኢባ > ο) ኮን ት፦ 


TABLE III 


DESIGN CONSTRAINTS 


Constraint 


Static Strength at Midspan 
Static Strength at x = 0 
Deflection 

Buckling due to Torque 
Buckling due to Compression 
Column Buckling 

Dynamic Strength (Fatigue) 
Speed 
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ЕХАЯРГЕ 1. : INETIAL DESIGN 
SHAFT DIMENSIONS : 
RI = 0.10000E Οἱ TH = 9.10000E Οἱ SL = 0.120008 
DESIGN CONDITIONS : 
ECCENTRICITY = 0.10000: 00 
END CONDITIONS 2 PINNED-PINNED 


LOADS + 

HORS EPOWER RPM PRESSURE AXIAL LOAD 
0.15000E 03 0.30000Е 03 0.0 0.0 

ANALYSIS RESULTS 

TORQUE(1) 0.31509E 05 
VOLUME 0.1131ЈЕ 04 
WEIGHT 0.31893F 03 
MOMENT ጋዮ INERTIA 0.11781E 02 
POLAR MOMENT JF INERTIA 0.23562E 02 
CRITICAL- SPEZI 0.14828E 04 
CRITICAL AXTAL LOAD 0.1875%Е 06 
CRITICAL BUCXLING TORQJE 0.47%5ЭЕ 08 
CRITICAL BUCKLING STRESS 0.905%9Е 07 
DEFLECTION 0.2658 +Е-01 
ВЕМОІМС MOMENT 0.50477Е 04 
SHEAR FORCE 0.2071%Е 03 
AXIAL STRESS AT A 0.85692E 03 
AXIAL STRESS AT B -0.85692€ 03 
HOOP STRESS 0.25013E 01 
TORSTINAL SHEAR STRESS 0.26738E 0% 
TOTAL SHEAR STRESS 0.27178Е 0% 
MEAN AXIAL STRESS 0.%%*T59E 02 
ALTERNATING AXIAL STRESS 0.81216Е 03 
VON MISES AT MIDSPAN 0.%7095Е 04 
VON MISES AT X=0. 0.%7073Е 04 
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EXAMPLE lo. : OPTIMUM DESIGN 
SHAFT DIMENSIONS : 
RI: 0.25847: Ol TH s 0.lil24€ 00 SL = 0.12000E 
DESIGN CONDITIONS : 
ECCENTRICITY = 0.10000E 00 
END CINDITIONS = PINNED-PINNED 


LOADS : 

HORSEPOWER RPM PRESSURE AXIAL LOAD 
0.159008 03 0.300008 03 0.0 0.0 

ANALYSIS RESULTS 

TORQJE(1) 0.31503Е 05 
VOLUME 0.22145£ 03 
WEIGHT 0.62451E 02 
MOMENT OF INERTIA 0.64353E 01 
POLAR MOMENT JF INERTIA 0.12871Е 02 
CRITICAL SPEED O.24757F 04 
CRITICAL AXIAL LOAD 0.45137E 05 
CRITICAL BUCKLING TORQJE 0.31509E 06 
CRITICAL BUCKLING STRESS 0.74729E 06 
DEFLECTION 0.92 753E-02 
BENDING MOMENT 0.95479E 03 
SHEAR FORCE 0.39681E 02 
AXIAL STRESS AT ለ 0.39998E 03 
AXIAL STRESS AT B -0.39998< 03 
HOOP STRESS 0.51655Е 01 
TORSIINAL SHEAR STRESS 0.65980F 04 
TOTAL SHEAR STRESS 0.66%*12E 04 
MEAN AXIAL STRESS 0.75507E O1 
ALTERNATING AXIAL STRESS 0.39243E 03 
VON MISES AT MIDSPAN 0.11435Е 05 
УОМ МІ5Е5 АТ Х-0. 0.11503E 05 
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EXAMPLE 2. 2 INITIAL DESIGN 

SHAFT DIMENSIINS : 

RI = 0.10000: 01 TH = 9.10900E 01 SL = 0.12000€ 
DESIGN CONDITIONS : 

ECCENTRICITY = 0.10000E 00 

END CINDITIONS : PINNED-PINNED 

ርዐልዐ5 ፡ 
HORS EP OWER RPM PRESSURE AXIAL LOAD 
0.15000Е 03 0.30200Е 93 0.100008 34 0.0 


ANALYSTS RESULTS 


TORQJE (1) 

VOLUME 

WEIGHT 

MOMENT JF INERTIA 

POLAR MOMENT JF INERTIA 
CRITICAL SPE=D 

CRITICAL AXIAL LOAD 
CRITICAL BUCXLING TORQUE 
CRITICAL BUCXLING STRESS 
DEFLECTION 

BENDING MOMENT 

SHEAR FDRCE 

AXIAL' STRESS AT A 

AXIAL STRESS AT B 

HOOP STRESS 

TORSIINAL SHZAR STRESS 
TOTAL SHEAR STRESS 

MEAN AXIAL STRESS 
ALTERNATING AXIAL STRESS 
VON MISES AT MIDSPAN 

VON MISES AT X=0. 
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0.31500E 05 
O.1131)E 04 
0.31893E 03 
0.11781£ 02 
0.23562Е 02 
0.14828E 0 
0.18764E 06 
0.47459E 08 
0.90549E 07 
0.26584£-01 
0.504775 04 
0.2071%Е 03 
0.1855ЭЕ 0% 
0.14303E 03 
0.66917E 03 
0.26733E 04 
0.27178E 04 
0.104*+8E 04 
0.81216Е 03 
0.49093€ 94 
0.47545E 04 





EXAMPLE 2. :*  JPTIMUM DESIGN 

SHAFT OIMENSIONS : 

RI = 0.18625Е 01 TH = 0.25560E 00 SL = 0.120008 
DESIGN CONDITIONS : 

ECCENTRICITY = 0.10000Е 20 

ENO CONDITIONS : PINNEO-PINNED 

LOADS : 
HORS EPOWER RPM PRESSURE AXIAL LOAD 
0.15200Е 03 0.30000Е 03 90.10000Е 04 0.0 


ANALYSIS RESULTS 


TORQUE (1) 

VOLUME 

WEIGHT 

MOMENT JF INERTIA 

POLAR MOMENT OF INERTIA 
CRITICAL SPE?) 

CRITICAL AXIAL LOAD 
CRITICAL BUC<LING TORQUE 
CRIFICAL BUC<LING STRESS 
ШЕББЕСТІПМ 

BENDING MOMENT 

SHEAR FIRCE 

AXTAL STRESS АТ А 

AXTAL STRESS AT B 

HOOP STRESS 

TORSIONAL SHEAR STRESS 
TOTAL SHEAR STRESS 

MEAN AXTAL STRESS 
ALTERNATING AXIAL STRESS 
VON MISES AT MIDSPAN 

VON MISES AT XzO. 
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0.31500ET 05 
0.39958E 03 
0.11268Е 03 
0.555785 Οἱ 
0.13315E 02 
0.187542 04 
0.75144E 95 
0.23553E 07 
0.22602E 07 
0.16358£-01 
0.17475E O^ 
0. 72247E 02 
0.75T03€ 04 
0.64537E 04 
0.6548TE O^ 
0.503%%Е 04 
0.507738 04 
0. 70305E 04 
0.5402TE 03 
0.1125%Е 05 
0.10965Е 05 
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EXAMPLE 3. : 
SHAFT DIMENSIONS : 
RI = 0.10000% 01 TH = 
DESIGN CONDITIONS : 


ECCENTRICITY = G.10000E JO 
END CINDITIONS ፥ PINNED-PINNED 


LOADS +: 
HORSEPOWER RPM 


0.15200Е 03 0.30000E 03 


ANALYSIS RESULTS 


TORQUE (1) 

VOLUME 

WE IGHT 

MOMENT ጋዮ INERTIA 

POLAR MOMENT IF INERTIA 
CRITICAL SPEz:) 

CRITICAL AXIAL LOAO 
CRITICAL BUCKLING TORQUE 
CRITICAL BUCXLING STRESS 
DEFLECTION 

BENDING MIMENT 

SMEAR FORCE 

AXIAL STRESS АТ ለ 

AXIAL STRESS AT B 

BOOP STRESS 

TORSIINAL SHEAR STRESS 
TOTAL SHEAR STRESS 

MEAN AXIAL STRESS 
ALTERNATING AXIAL STRESS 
VON MISES AT MIOSPAN 

VON MISES AT Xz0. 
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INIFIAL DESISN 


0.10000E Οἱ SL = 0.12000E 


PRESSURE AXIAL LOAD 


0.20000E 04 


0.315026 05 
0.11310Е 04 
0. 31893= 03 
0.11781€ 02 
0.23562* 02 
0.14757Е 04 
0.18754E 06 
0.47469£ 08 
0.90549£ 07 
9.2681%Е-01 
0.541442 04 
0.2072)E 03 
0.7059Т7Е 03 
-0.11314E 04 
0.25913F 01 
0.26738Е 04 
0.27178Е 04 
-0.10512Е 03 
0.81216Е 03 
0.%75ТТЕ 04 
0.%7043Е 04 





EXAMPLE 3. : JPTIMUM DESIGN 

SHAFT DIMENSTINS : 

ШЕТ ӘУ ОТЕ οι ο ο 11119E 00 SL = 0.12000€ 
DESIGN CONDITIONS : 

ECCENTRICITY = 0.10000E 09 

END CONDITIONS : PINNEO-PINNED 
LOADS : 
HORSEPOWER RPM PRESSURE 
0.15200Е 03 0.30000Е 03 0.0 


AXIAL LOAD 
0.200008 O^ 


ANALYSTS RESULTS 


TORQUE (1) 

VOLUME 

WEIGHT 

MOMENT OF INERTIA 

POLAR MOMENT JF INERTIA 
ERTTIZCAL SPERO 

CRITICAL AXIAL LOAD 
CRITICAL BUCXLING TORQJE 
CRITICAL BUCCLING STRESS 
DERLE-TTON 

BENDING MOMENT 

SHEAR FORCE 

AXIAL STRESS АТ А 

AXIAL STRESS AT B 

HOOP STRESS 

TORSTINAL SHEAR STRESS 
TOTAL SHEAR STRESS 

MEAN AXIAL STRESS 
ALTERNATING AXIAL STRESS 
VON MISES AT MIDSPAN 

VON MISES AT X=0. 
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0.31500E 05 
0„22185=Е 03 
0.6256%Е 02 
0.64766E 01 
0.12953Е 02 
0.245357 04 
0.45241Е 05 
0.31500£ 06 
0.745275 06 
0.93 752Е- 02 
0.98171Е 03 
0.39758Е 02 


20.67221E 03 


0525218355 01 
0.65705Е 04 
0.66136Е 04 


-0.1063ТЕ 04 


0.39151Е 03 
0.114732 05 
0.11455E 05 
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EXAMPLE 6. + INITIAL DESIGN 


SHAFT DIMENSIINS : 


RI = 0.100002 01 TH = 9.10009E Ol SL 


DESIGN CONDITIONS ፡ 
ECCENTRICITY * 0.100008 00 


END CINDITIONS < PINNED-PINNED 


LOADS : 


HORSE? IWER RPM PRESSURE 


0.15090Е 03 0.30000Е 0% 0.0 


ANALYSIS RESULTS 
TORQUE (1) 
VOLUME 
WEISAT 
MOMENT OF INERTIA 
POLAR MOMENT JF INERTIA 
ወው! ፐ ላር 5РЕЗЭ 
CRITICAL AXIAL LOAD 
CRITICAL BUCKLING TORQUE 
CRITICAL BUCX<LING STRESS 
DESLE-TIDN 
BENDING MOMENT 
SHEAR FORCE 
AXIAL STRESS AT A 
AXIAL STRESS AT B 
НООР STRESS 
TORSIINAL SHEAR STRESS 
TOTAL SHEAR STRESS 
MEAN AXIAL STRESS 
ALTERNATING AXIAL STRESS 
VON MISES AT MIDSPAN 
VON MISES AT X=0. 


43 


0.315008 
0.11310E 
0.31893E 
0.11781€ 
0.23562Е 
0.14828E 
0.1876%Е 
0.47459€ 
0.90549E 
0.24000E 
0.23801E 
0.627225 
0.40405E 
-0.40%06Е 
0.25013E 
0.267385 
0.13337€ 
0.40397E 
0.81200Е 
0.404945 
0.23100€ 


= 0.12000E 


AXTAL LIAD 
0.0 


04 
04 
03 
02 
02 
04 
06 
08 
07 
02 
08 
06 
07 
07 
03 
03 
06 
07 
03 
07 
06 





EXAMPLE 42 2 IPTIMUM DESIGN 
SHAFT DIMENSIINS > 
αν ορ 01 Tri s 9236963E-OL SU = 0.120006 
DESIGN CONDITIONS : 
ECCENTRICITY = 0.100008 00 
END CONDITIONS : PINNED-PINNED 


LOADS : 
HORSEPOWER RPM PRESSURE AXIAL LOAD 
0.15000E 03 0.30000Е 04 0.0 0.0 
ANALYSIS RESULTS 
TORQUE(1) 0.31500E 04 
VOLUME 0.17831E 03 
WEIGHT 0.502835 02 
MOMENT OF INERTIA 0.30%15Е 02 
POLAR MOMENT IF INERTIA 0.60832E 02 
CRIFIZAL SPE:) 0.60005Е 06 
ο απ ο ΑΧ LOAD 0.39475£ 05 
CRITICAL BUC<LING TORQJE 0.31%9%Е 05 
CRITICAL BUC<LING STRESS 0.10432£ 06 
DELETION 0.43922E-01 
BENDING MOMENT 0.76214E 04 
SHEAR FORCE 0.84TT3E 03 
AXIAL STRESS AT A 0.16079E 04 
AXIAL STRESS AT B -0.16079E 04 
HOOP STRESS 0.29628E 04 
TORSIINAL SHEAR STRESS 0.33228£ 03 
TOTAL SHEAR STRESS 0.14733Е 04 
MEAN AXIAL STRESS 0.14488E£ 04 
ALTERNATING AXIAL STRESS 0.15913E 03 
VON MISES AT MIDSPAN 0.2632 TE 04 
VON MISES AT X=0. 0.391032 04 
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EXAMPLE 5. : INITIAL DESIGN 
SHAFT DIMENSI 3 NS : 
RIs 0.10000 Οἱ ТН = 0.10000Е 01 οἱ 
DESIGN CONOITI3NS : 
ECCENTRICITY = 0.10000= 00 
END CONDITIONS : PINNED-PINNEO 


LOADS : 
HORSEPOWER 
0.15000E 03 


RPM PRESSURE 
0.30000E 04 0.100002 04 


ANALYSIS RESULTS 


TORQUE (1) 0.31509€ 
ከ ነገር 0.11312Е 
WEIGHT 0.31893€ 
MOMENT OF INERTIA 0.117815 
POLAR MOMENT JF INERTIA 0.23552€ 
GRITIZAL SPEED 0.14828Е 
CRITICAL AXIAL LOAD 0.1875*E 
CRITICAL BUC<LING TORQJE 0.%7%59Е 
CRITICAL BUC<LING STRESS 0.905%9Е 
DEFLECTION 0.24000E 
8ENDING MOMENT 0.23801E 
SHEAR FORCE 0.62722Е 
AXIAL STRESS AT А 0.40416Е 
AXIAL STRESS AT 8 -0.40396Е 
HOOP STRESS 0.91680E 
TORSIINAL SHEAR STRESS 0.26738E 
TOTAL SHEAR STRESS 0.13334E 
MEAN AXIAL STRESS 0.49407E 
ALTERNATING AXIAL STRESS 0.81209€ 
VON MISES AT MIDSPAN 0.404112 
VON MISES AT X=9. 0.231098 
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= 0.12000E 


AXTAL LOAD 
0.0 


04 
04 
03 
02 
02 
04 
06 
08 
07 
02 
08 
06 
07 
07 
03 
03 
06 
07 
03 
07 
06 





EXAMPLE 5. : OPTIMUM DESIGN 

SHAFT DIMENSI INS 3 

፳፲ = 0.608442 01 ТН = 0.51266Е 00 SL = 0.120005 
DESIGN CONDITIONS ፡ 

ECCENTRICITY = 0.10000E 00 

END CINDITIONS 2 PINNED-PINNED 
LOADS : 
HORS EPOWER RPM PRESSURE AXTAL LOAD 
0.15200Е 03 0.309000Е 0% 0.100005 04 0.0 


ANALYSIS RESULTS 


TORQUE(1) 0.31509Е 04 
VOLJME 0.29521E 04 
WEIGHT 0.832598 03 
MOMENT OF INERTIA 0.50353E 03 
POLAR MOMENT JF INERTIA 0.10071E 04 
CRITICAL SPEED 0.60003E 04 
CRITICAL AXIAL LOAD 0.6535ТЕ 06 
CRITICAL BUCCLING TORQJE 0.34401£ 08 
CRITICAL BUCKLING STRESS 0.16557E 07? 
DEFLECTION 0.43925E-01 
BENDING MOMENT 0.12619Е 06 
SHEAR FORCE 0.14035E 05 
AXIAL STRESS AT A 0.11609E 05 
AXIAL STRESS AT B 0.82525Е 04 
HOOPRESTRES S 0.12589E 05 
TORSIINAL SHEAR STRESS 0.209%85 02 
TOTAL SHEAR STRESS 0.11620Е 04 
MEAN AXTAL STRESS O.11443E 05 
ALTERNATING AXTAL STRESS 0.16609E 03 
VON MISES AT MIDSPAN 0.12129E 05 
VON MISES AT X=0. 0.12749E 05 
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EXAMPLE 6.2 : INITIAL DESIGN 


SHAFT DIMENSIINS : 


RI = 0.10000: 01 TH = 9.10000E OL SL 


DESIGN CONDITIONS : 
ECCENTRICITY = 0.10000Е 00 


END CONDITIONS < PINNED-PINNED 


LOADS : 


HORSEPOWER RPM PRESSURZ 


0.15000E 93 0.30200Е 0% 0.0 


ANALYSIS RESULTS 
TORQJE (1) 
VOLUME 
WEIGHT 
MOMENT JF INERTIA 
POLAR MOMENT JF INERTIA 
"የ cal ου -ὂ 
CRITICAL AXIAL LOAD 
CRITICAL BUCKLING TORQUE 
CRITICAL BUCKLING STRESS 
DEFLECTION 
BENDING MOMENT 
SHEAR FORCE 
AXIAL STRESS AT A 
AXIAL STRESS AT B 
HOOP STRESS 
TORSTINAL SHEAR STRESS 
TOTAL SHEAR STRESS 
MEAN AXIAL STRESS 
ALTERNATING AXIAL STRESS 
VON MI SES AT MIDSPAN 
VON MISES AT X=0- 
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0.315002 
ο 0515 
0.31893€ 
0.11731€ 
0.23552€ 
0.14757E 
0.13754E 
0.47%5ЭЕ 
0.905%9Е 
0.24000E 
0.56420€ 
0.62722Е 
0.95779Е 
-0.95733E 
0.25013E 
0.267336 
@.1333TE 
0.95771Е 
0.812508 
0.95785€ 
0.23100E 


= 0.12090Е 


AXIAL LOAD 
0.20000E 04 


04 
04 
03 
02 
02 
04 
06 
08 
07 
02 
08 
06 
07 
07 
03 
03 
06 
07 
03 
07 
06 





EXAMPLE 6. : 
SHAFT DIMENSTINS : 


JPTIMUM DESIGN 


RI = 0.63959E O1 TH = 9.36942E-01 SL = 0.12000E 


DESIGN CONDITIONS : 


ECCENTRICITY = 0.150000Е 90 
END CINDITIONS : PINNED-PINNED 


LOADS : 
HORSEPOWER RPM 


0.15200Е 03 0.30000Е 04 


ANALYSIS RESULTS 


TORQJE(1) 

VOLUME 

WEIGHT 

MOMENT ЭР INERTIA 

POLAR MOMENT JF INERTIA 
CoPII AL 5ΡΕΞη 

CRITICAL AXIAL LOAD 
CRITICAL BUC<LTNG TORQUE 
"ከ ከ! ለእ. BUC<LING STRESS 
DEFLECTION 

BENDING MOMENT 

SHEAR FIRCE 

AXIAL STRESS AT A 

AXIAL STRESS AT B 

НООР STRESS 

TORSTINAL SHEAR STRESS 
TOTAL SHEAR STRESS 

MEAN AXIAL STRESS 
ALTERNATING AXIAL STRESS 
VON MISES AT MIDSPAN 

VON MISES AT X=0. 
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PRESSURE AXIAL LOAD 


0.20000E 04 


0.31500E 04 
0.17856= 03 
0.50382F 02 
0.396292 02 
0.61258Е 02 
0.60051Е 04 
0.3955ТЕ 05 
0.31490E 05 
0.10400E 06 
0.43813E-0l 
0.17025E 05 
0.84895£ 03 
0.22325E 04 
-0.49193E 04 
0.29778Е 04 
98.2330T93E 03 
0.14712E 04 
0.209733Ε 04 
0.15872Е 03 
0.69313E 04 
0.39193E 04 
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ERIMPLE T. 2 INITTIA 
SHAFT DIMENSTONS : 
КІ = 0.10000= οἱ TH = 0.1000 
DESIGN CONDITIONS : 
ECCENTRICITY = 0.10000€ 90 
END CINDITIONS : PINNED-PINNED 
LOADS 5 


C DESIGNY 


SP OI SL 


HORSEPOWER RPM PRESSURE 


0.15200Е 03 0.30090€ 03 9.0 
0.15300Е 03 9.30000Е 0% 59.0 


ANALYSIS RESULTS 
TORQUE(1) 
TORIJE (2) 
VOLUME 
WEIGHT 
MOMENT IF INEITIA 
POLAR MOMENT IF INERTIA 
CRITICAL 5ΡΕΞΠ 
CRITICAL AXIAL LOAD 
CRITICAL BUCKLING TORQUE 
CRITICAL BUCKLING STRESS 
DEFLECTION 
BENDING MOMENT 
SHEAR FIRCE 
AXIAL STRESS AT A 
AXIAL STRESS AT B 
HOOP STRESS 
TORSTONAL SHEAR STRESS 
TOTAL: SHEAR STRESS 
MEAN AXTAL STRESS 
ALTERNATING AXTAL STRESS 
VON MISES AT MTDSPAN 
VON MISES AT X=0. 
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0.31500E 
0.31500E 
O SIDE 
0.31893€ 
0.117381: 
0.23562E 
O.1476TE 
0.18754E 
0.47469E 
0.90543ЭЕ 
0.2490IE 
0.56%20Е 
0.627225 
0.95773€ 
-0.95783E 
0.25013E 
0.26733€ 
0.133378 
0.957{1ЄЕ 
0.81259€ 
0.95785Е 
0.23100Е 


- 0.12000Е 


AXTAL LOAD 
0.0 
0.20000Е 04 


05 
04 
04 
03 
02 
02 
04 
06 
08 
07 
02 
08 
06 
07 
07 
03 
03 
06 
07 
03 
07 
06 





ΠΕ .: 
SHAFT OIMENSIINS : 
АІ = 0236355215 01 TH = 
DESIGN CONDITIINS : 
ECCENTRICITY » 0.100007 


END CONDITIONS : PINNED-PINNED 


LOADS : 
HORSEPOWER RPM 


0.15000Е 03 0.30000E 03 
0.152005 03 92.30000Е 04 


ANALYSIS RESULTS 


TORQJE(1) 

TORQJE (2) 

VOLUME 

WEIGHT 

MOMENT IF INRTIA 

POLAR MOMENT IF INERTIA 
CRITICAL SPEED 

CRITICAL AXIAL LOAD 
CRITICAL BUC<CLING TORQUE 
CRITICAL BUCKLING STRESS 
DEFLECTION 

BENDING MOMENT 

SHEAR FIRCE 

AXIAL STRESS AT A 

AXIAL STRESS AT B 

HOOP STRESS 

TORSTINAL SHZAR STRESS 
TOTAL SHEAR STRESS 

MEAN AXIAL STRESS 
ALTERNATING AXIAL STRESS 
VON MISES AT MTDSPAN 

VON MISES АТ Х=0. 


S 


JPTIMUM DESTSN 


0.93009Е-01 SL = 0.12000€ 


PRESSURE AXIAL LOAD 


0.0 
0.20000E 0% 


0.315998 05 
0.315008 04 
0.44893E 03 
0.12660Е 03 
0.76660E 02 
05153325 03 
0.60002Е 04 
0.99390E 05 
0.31571Е 06 
0.26122E 06 
0.43932E-01 
0.42893E 05 
213546 05 
0.30733E 04 

-0.41425E 04 
0.29828E 04 
DNI EC US 
0.12736E 04 
0.29136Е 06 
ር оз 
0.620215 04 
0.37098E 04 





C BI X ን ӨЗНӘНІ SNI aU 30915ክ1 


ga“? gua ga'g Bh ° B “ዳይ. gua 





umurr 3do 


"ur "πο 00S 


uorboNy әтатвеәд 


σας 


/. ጋገ)ዘ8።ጋ 


ፍኤ3ዘፈን | ዘ 


( BX )S3DNI 
54 





B. OBSERVATIONS AND CONCLUSIONS 

With all other parameters held constant, an increase in 
speed drives the design to a very large inside radius while 
reducing the thickness considerably. This trend is so 
strong that a reasonable lower bound on the thickness must 
be given to prevent numerical ill-conditioning. The intro- 
duction of eccentricity increases the effect of rotation on 
the stresses and the optimization produces slightly larger 
values for the design variables. This effect is reduced, 
however, by the constraint on the deflection. 

As the internal pressure is increased, the tendency is 
to produce an optimized radius that is smaller, and a larger 
thickness. This may be explained by the direct effect of 
the internal pressure on the hoop stress which in turn would 
require a larger thickness to support the higher stress. 

The effect of axial load is reflected on the deflection. 
Again, if very small deflections are allowed, the axial load 
effects are reduced. The major parameter that is affected 
is the axial stress. 

A notable characteristic shown by the two variable 
function space diagrams is the tendency of the strength 
constraints to curl back as the radius is increased. This 
signifies the existence of an upper limit on size for a given 
design load. The increase in mass, especially at high speeds, 


increases the stresses considerably leading to violation of 


25 





the strength constraints. This phenomenon is primarily a 
result of the eccentricity of the mass but also is due to 
the static weight of the shaft itself. 

It is also noted that these strength constraints generally 
run close to each other in an almost parallel manner. The 
dynamic or fatigue constraint maintains a position such that 
coming from the feasible range towards the infeasible area, 
it is activated before the static constraints are encountered. 
In effect, the latter constraints are redundant in rotating 
shafts. However, these are kept in the analysis to allow 
solution of static problems. 

The versatility of the program should not be lost in the 
emphasis on drive shaft applications. With the appropriate 
loadings, one may use it to solve structural problems like 
columns, beams and pressure vessel designs. The objective 
function need not be confined to the volume. For example, 
given a specific initial dimensions, one may minimize the 
Young's modulus to determine whether or not a certain 
material may be utilized; or to determine how low the modulus 
could get before the shaft fails. If the outer dimension 
is critical to the design, one may set a specific value on 
the diameter and let all other parameters vary during 
optimization. 

In summary, numerical optimization provides an efficient 


and effective way of improving shaft designs. The analysis 
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used in this exercise is a first step when compared to the 


boundless refinements in techniques that are possible. 
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V. RECOMMENDATIONS: FOR FUTURE INVESTIGATIONS 


The study has shown the feasibility of using numerical 
optimization techniques in the design of drive shafts within 
the limitations imposed in the analysis. Further studies on 
the same design field may be pursued by eliminating some of 
these limitations. For example, the unfeasible region defined 
by the speed constraint extends only to some discrete distance 
and a supercritical feasible design area does exist. It will 
be interesting to develop optimization techniques to accommo- 
date this disjoint design field and to investigate the pecu- 
liarities this will introduce to the analysis. 

Another investigation may be concerned with shafts with 
variable thickness and radius along its length. The inclusion 
of radial stresses in the calculations, the use of finite 
element methods in three dimensional stress analysis of the 
problem, and vibration considerations are but some of the 
aspects of research that may be useful in future analyses of 
composite materials. 

The challenge of including the unique material character- 
istics of composites in the design will invariably stimulate 
a lot of young minds. However, the handling of these 
direction oriented properties in stress analysis have already 


been investigated [6 through 12], and the tools required to 
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carry out the program are available. COMAND [13] offers a 
simplified composite analysis procedure for panels. Such 
procedures, adapted with the program used in this investi- 
gation could result in a simplified approach to composite 
drive shaft analysis. Іп this procedure, failure of the 
ΤτΕ 15 assumed if a single ply; fails and the entire 
analysis is performed in a ply by ply basis. Since the 
stress elements have to be in a differential level, the 
composite material is assumed to be homogenous and uni- 
directionally isotropic at the differential range. It will 
also be necessary to assume that stresses throughout the 
thickness of the ply is uniform and equal to the stresses 
at the outer element of the ply. Analysis may then be 
performed in the classical procedure, solving for the 
stresses along the major shaft axes. The material properties 
to be used in the calculation must be equivalent values 
derived from the directional properties using transformation 
of axes equations. References 10 and 13 present procedures 
to calculate the corresponding material properties, stresses 
and strains along the major axes of the shaft and the fiber 
orientation coordinates. 

More involved investigations could deal with failure 
modes of the whole shaft instead of a single ply, three 
dimensional analysis of plys including transverse shear 


stresses, delamination problems, thermal effects, and others. 
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The field of composites also opens up another group of 
considerations like the effects of stacking different ply 
orientations at different sequences, unsymmetrical laminates, 


unequal thicknesses of plys and other factors relevant only 


to composites. 
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APPENDIX A 


The FORTRAN program used in this investigation is 
described here. Subroutine ANALIZ provides the basic 
analysis used in the optimization. It reads the initial 
design description and calculates the values of the objec- 
tive function, constraints, and all other parameters 
necessary to solve the analysis problem.  COPES-CONMIN 
updates the design to minimize/maximize the objective 
function, iterating until no further improvement in the 
objective function is possible without violating one of 
the constraints. This process is shown in the flow diagrams. 

The global catalog lists the location, FORTRAN name, 
mathematical symbol and description of the parameters used 
in the optimization. These parameters are contained in the 
labeled COMMOM block, GLOBCM. 

There are at least five input cards to describe the 
initial design. The first card contains the initial 
dimensions--inside radius, thickness and shaft length. The 
second card describes the material used. Young's modulus, 
shear modulus, Poisson's ratio, weight per unit volume, 
yield strength, ultimate strength and compressive strength 
are read here. The endurance limit is internally calculated 


using the endurance limit factors read from the third card. 
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The size factor, CB, is internally updated as the shaft 
diameter is varied during optimization. 

The fourth and fifth cards specify the design and load 
conditions. The parameter, LC, allows the user to choose 
the end conditions applicable to his requirements. Pinned- 
pinned and clamped-clamped analyses are available while 
clamped-pinned and cantilever end conditions are yet to be 
developed. The parameter LCC allows multiple loading design. 
One can thereby design a shaft that is capable of supporting 
LCC different load combinations. 

The magnitude of mass imbalance may be read as a specified 
distance from the center of the shaft, EO, or as a specified 
fraction, EC, of the shaft radius. The safety factor 
generally used in the analysis is also specified in the fourth 
card. 

The fifth card specifies the horsepower, shaft speed in 
RPM, internal pressure and axial load for each loading. 
Additional cards are used up to LCC number of loadings. The 
torque is calculated internally in inch-pounds. The follow- 
ing table summarizes the required information for subroutine 


ANALIZ. 
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TABLE IV 


INPUT CARDS 





CARD INFORMATION 

1 RI, TH, SC 

2 E, G, PR, SPWT, I/S 
SU, YC 

3 СИ ен, СЕК р СЕ, 
СЕ 

4 UG eee, κο EC, SF 

5 HP, RPM, P, F 


FORMAT 
δει ο. 2 
δι 0.2 
SF10 2 


21975500272 
ὄν 7 


Subroutine BEND calculates the bending moment, deflection 


ЕШ ምጋ ከ r force at the critical cross-sections of the beam. 
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SUBROUTINE ANALIZ (ICALC) 
FLOW DIAGRAM 


ICALC = 1 


READ INITIAL DIMENSIONS 
READ MATERIAL PROPERTIES 
READ DESIGN CONDITIONS 
PRINT ALC TNPUT 


ICALC = 2 


CALCULATE GEOMETRIC PARAMETERS 
CALCULATE OBJECTIVE FUNCTION 
CALL SUBROUTINE BEND 
. CALCULATE MAXIMUM DEFLECTION 
CALCULATE BENDING MOMENTS 
CALCULATE SHEAR FORCE 


CALCULATE STRESSES AT CRITICAL STRESS 
ELEMENTS 


CALCULATE CONSTRAINTS 


ICALC = 3 


PRINT ANALYSIS RESULTS 


64 





OPTIMIZATION FLOW DIAGRAM 


TE 
SUBROUTINE SUBROUTINE 
ANALIZ gen BEND 
COPES CONMIN 
...::2=፦-= ሙ=: 

OUTPUT 
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Global 
Location 


Ww N FP 


5-14 
15-24 
25-34 
35-44 
45 
46 
47 
48 


49 


50 
51 
52 


53 
54 


Бо 
56 
27 
58-67 
68-77 


78 
79 


TABLE V 


GLOBAL CATALOG OF PARAMETERS 


FORTRAN 
Name ` 
ВІ 

TH 

SL 

VOL 

το (10) 
P(10) 

F CEO) 
1015 πο 
E 

G 

PR 

NS 


1C 


SU 
RO 
SPWT 


AI 
AJ 


HS 

AST 
ASTN 
ASC (10) 
SALT (10) 


SMEAN 
SST 


Math 
Symbol 


[5- 


< GQ tj € hj 'U Hd s ኤፎ rr H 


ሀ 
= 


un 
Q 
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Definition 


Inside Radius 
Thickness 

Shaft Length 
Volume 

Torque 

Internal Pressure 
Axial Load 
Shaft Speed 
Young's Modulus 
Shear Modulus 
Poisson's Ratio 


Yield Strength 
(Tension) 


Yield Strength 
(Compression) 


Ultimate Strength 
Outside Radius 


Weight Per Unit 
Volume 


Moment of Inertia 


Polar Moment of 
Inertia 


Deep 5ῬΕΤΕΞΞ 

Axial Stress at A 
Minimum Axial Stress 
Axial Stress at B 


Alternating Axial 
Stress 


Mean Axial Stress 


Torsional Shear 
Stress 
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8l 
82 


83 


8-93 
94-103 
104-113 
114 

12155 

1 PG 

11.7 
118 


119 
ЕЛ 01-129 
ሪጋ 


SSV 


SSMX 
ECRI 


SCSX 


VMSI (10) 
VMS2 (10) 
VMS3 (10) 
WM 

BM 

BMN 

V 

ЕС 


SF 
DEL (10) 
GC (10,10) 


Q 


Q Q 
O 


mid 


= E 


max 


= 


min 


Max 


wire 
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Shear Stress due to 
Sheer Force 


Total Shear Stress 


Critical Torque for 
Buckling 


Critical Stress in 
BUCkling 


Von Mises (mean) 

Von Mises at x= 0 

Von Mises at Midspan 
Critical Speed 

Maximum Bending Moment 
Minimum Bending Moment 
Shear Force 


Critical Column Buckling 
Load 


Safety Factor 
Deflection 


Constraints 
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APPENDIX B 


The derivation of the equations for the deflection, 


bending moments and shear force at the critical cross- 


. " | | . ም | 
Ectrons cof the shaft are presented here. Pinnedrpinned 


and fixed-fixed end conditions are analyzed. 


1. Pinned-Pinned 





Let y = do sin => 
ао = ам (у+е) πό 
dQ = ШІ (y+e) dx 
Let К, = эрине 
dQ = к, ια. sin - + e] dx 
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CE) 


E 2) 


(13) 





With the effect of weight: 


dQ = [K, (d sin Шалы 1 መመር 


7 2] dx (1.4) 


where K, = pAg = numerically equal to weight per unit length 


| ደ 2 3,፻ 
Q = | dQ = Kil zer] ES Q (125) 
Moment at x = Xo 
Da O (p= x) + CIE S $93 
хо 2 О О : 
О 
а а 
_ о Lc MC: OTE _О 
МЕ = + Ky ЖІ + 2) + [5 Ку? (7 + е) |х (1.7) 
к.е 
1 2 2 TX 
τ». - K (5) Sin Бои 
m 2 
2 2 


This represents the moment at any value of x. If the axial 


load, F, is included, the total bending moment is 


MTOTAL ^ Hx. t Mp = Ely 


where M. = Fy = F (d, sin =) (1.8) 


Integrating twice: 


K.e 3 
zanl y 3 4 Bar. TX 
Εν = 547 [2 x 22x + x ] Б Ky do . Sin 2 (1.9) 
2,2 E. νι de 1 
+ P d. (—) San ШЕ о [2 EX ] 
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5 K ат" Beker ae 


2,4 DONNEES Ὁ 

ን = do m ο u 384 384 

a - 50% o 10) 
Ο 384 572 i 


LSA 
Е = K = s =) 


This is the deflection when the weight effect is subtractive. 
The maximum deflection occurs when the weight tends to increase 


the deflection. 


" _ 554 | куе + к. 
тах 384 2 2 05-22 
EI - K, (7) _ F (=) 


] E 


Substituting (1.10) in (1.7) and (1.8), the bending moment 
at x = P is 


K,,2 


8 


= E 2,2 = 
M in ፦ K; ls + do С) ] + Fd 





(Laa) 


Again the maximum moment occurs when the weight effect is 








additive 
2172 
E e 2 2 
— Κι ο 4 que (=) ] + Fd + 5 (1013) 
From Equation (1.5), the maximum shear force VMAX is equal 
Q 
to 2 Or 
g К.е? Ku 
1 "max 1 2 
VMAX = s: 2 sl: SON (1.14) 
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2. Clamped-Clamped 


A similar derivation procedure for clamped-clamped ends 


using a shape function: 





95 ር በበረታ 
y= => [1 - cos 7 ] 
and | 
do do ATX 
dQ = K, [= - - cos — + e] ах 
do 
О = dQ = Ki? [I> + e] 
G 
Moment at x = Xo 
X 
ዕ О 
M= M -=x + | (x - x) dQ 
О 2 
О 
K,2 do Ж X d 
Ὁ... η — 
О 
ха ха 
О О 20 
+ xe => በጋ! a a xe 
M at any value of x: 
K $? а K 
ን ሸበ » AO ΝΗῚ 
Py > ML OM, =з ME E 
к.а 12 
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O 

2 COS 
md 

+ e) a 


(231) 


(2%2) 


(225) 


2тх 


(29) 





Integrating once: 


























K; do «3 x^12 куа 2х 
= zx M ο ας u A 
Ely = M, ž A ISS 5—) 5 
8T 
3 
τ ο. 
Т ТО эе pu 
16T 
y' = 0 ሕቲ ሄ= 0). P a c, = 0 
M =K. 2° [o . 9 у со, 
O ጊ ሪ4 r2 2 
8T 
Integrating once more: 
nee, 
Ку 2 24 12 24 
ae" ATX 
= cos —— + С 
3244 ie 2 
μα 
y = 0 at x = 0, C. = 
? 327" 
EI a κ 2 ae" ATX 
IL + e) (E) (A) + (1 - соз ---) (27-6) 
K 2 24 4 ዷ 
Т Su 
у= а _ at x = 6/2: 
а 2 72 545 
7, ደ 
| 16т 
4 4 
Elda, _ d ‚ е0“ f αμ 
Ky 768 384 167“ 


(5 





B oU —— — — omm: Ей) 


With axial load F: 

















= = E x o А 27x 
Moment due to F M. М. "Ys M 2 (1 cos i 
Integrating once: (2.8) 
i Fd 
= » О ሺ ፣ ከ 2x 
Ely, = М. x =~ (x 57 Sin mo кт 
Ур =0 at x= 0, C1 = 0 
: Fd, 
Y, -Üatx-e,M,- —— 
where М. = Moment at the ends due to the ахіа1 load 
ae x? > x? 12 απ. 
E e Τ᾿; ΤΊ: 
ፐ 
ға 44 
Yp = 0 at x = 0, с. = 5 
8T 
rà 1^ ATX 
Ур = 2 (1 - COS 2 ) (2.9) 
St EL 
At x = 1/2 
ға μ΄ 
vy. = 
ч 4ፐ“ጄ፲ 
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Eom Equation (2.6) 


2 2 














1 а 2 D x D 
πι (53) (x-2) + e (53) (፳ > ደ) 
Е E 27 
+ — (1 — cos 35] 
232m 
At x = PE 
4 4 
Fest a Ке v 
1 ዕ 1l O 
y. = =— |= + e] + 
м 384EI 2 1є6т^®т 
do E ሃጅ Ы Yg 
e куя“ кул“ ТАС 
= ( + m + ) ኋ + (210) 
mom 768EI 16т^ЕТ О 384EI 
4 4 
2 K. 24K,2 
er da = (292 ን ከ. 
О 2 4 ο T 
TT 2 T 
4 
Ғы |i (2.11) 
О 2 
QNM 24 96፻2 
384EI K; ዷ (5 + 74) - > 
T 
When the effect of weight is considered: 
4 4 
a куех £ КЁ (2.12) 
πο -. ш-н 
O  384EI - 0.74638K, 2° - 96(2) ፻ 
where K, = pAg, numerically equal to weight per unit length. 
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πο, 


11: 


12. 


їз; 
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